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We present a simple analytical model and an exact numerical study which explain the role of 
roughness on different length scales for the fluid contact angle on rough solid surfaces. We show 
that there is no simple relation between the distribution of surface slopes and the fluid contact angle. 
In particular, surfaces with the same distribution of slopes may exhibit very different contact angles 
depending on the range of length-scales over which the surfaces have roughness. 


Controlling surface superhydrophobicity is of utmost 
importance in a countless number of applications [ij. 
Exarmies are anti-icing coatings [3, , friction reduc¬ 

tion Q , antifoggiM properties [J, antirefiective coat¬ 
ings [7[, solar cells (8|, chemical microreactors and mi¬ 
crofluidic microchips (9|, self-cleaning paints and opti¬ 
cally transparent surfaces (lol - [l^ . Water droplets on 
superhydrophobic surfaces usually present very low con¬ 
tact angle hysteresis, i.e. very low rolling and sliding 
friction values, which make them able to move very eas¬ 
ily and quickly on the surface, capturing and removing 
contamination particles, e.g, dust. Surfaces with large 
contact angles, now referred to as exhibiting the “lotus 
effect”, are found in many biological systems, such as 
the Sacred Lotus leaves [l^ 141, water striders [l^, or 
mosquito eyes @ , where the presence of surface asperities 
cause the liquid rest on the top of surface summits, with 
air entrapped between the drop and the substrate. This 
type of “fakir-carpet” configuration is referred to as the 
Cassie-Baxter state [l^. However, depending on the sur¬ 
face chemical properties and its micro-geometry, a drop 
in a Cassie-Baxter state may become unstable when the 
liquid p ressure increases above a certain threshold value 


1IH25|. 


The increase of liquid pressure may occur during drop 
impacts at high velocities [2^. When this happens the 
liquid undergoes a transition to the Wenzel state 27j| . 
which makes the droplet rest in full contact with the 
substrate [2l|, [2^. This transition to the Wenzel state 
may be irreversible because of the energy barrier the drop 
should overcome to come back to the Cassie-Baxter state 
17 . 2^, 28, 2^. However, the presence of multi-scale or 


hierarchical micro-structures may favour the transition 
back to the Cassie-Baxter state |3Cll432i. o r the stabiliza¬ 
tion of the ‘fakir-carpet’ configuration [20j, thus explain¬ 
ing why many biological systems present such a multi¬ 
scale geometry M. 

Some studies |35l - l40l| , have shown that many randomly 
rough surfaces possess super water-repellent properties 
with contact angles up to 174°, which could explain why 
biological systems use such hierarchical structures to en¬ 
hance their hydrorepellent properties Only few theo¬ 


retical studies focus on this aspect of the proble ml35l . l36 . 
Ifll - li^ . In particular, it has been suggested 33-40| that 
for randomly rough surface the critical parameter which 
stabilizes the Wenzel or Cassie state is the so called Wen¬ 
zel roughness parameter r. Given Young’s contact angle 
0 the Cassie state is stable when r > — [cos0]~^, on the 
other hand when r < — [cos0]~ the low energy state is 
the Wenzel state. In this Letter we present exact numeri¬ 
cal results, and the first analytical theory, which show the 
fundamental role of roughness on many length scales for 
fractal-like surfaces, in generating large contact angles. 

Consider a fluid droplet on a perfectly flat substrate. 
If the droplet is so small that the influence of the gravity 
can be neglected, the droplet will form a spherical cup 
with the contact angle 9. In thermal equilibrium the 
Young’s equation is satisfied: 


7SV = 7SL + 7LVCOS0 


(I) 


Consider now the same fluid droplet on a nominal flat 
surface with surface roughness. If the wavelength Aq of 
the longest wavelength component of the roughness is 
much smaller than the radius of the contact region, the 
droplet will form a spherical cup with an (apparent) con¬ 
tact angle 9q with the substrate, which may be larger 
or smaller than 9 depending on the situation. In this 
case the contact angle 9o will again satisfy the Young’s 
equation, but with modified solid-liquid and solid-vapor 
interfacial energies 7 gL and qgy (see also Ref. 0): 

7sv = 7 sl + 7lvcos6»o (2) 

If the fluid makes complete contact at the droplet- 
substrate interface (Wenzel model), then 7 gy = rjsv 
and 7gL = r7sL, where r = Htot/^o is the ratio between 
the total surface area of the substrate surface, and the 
substrate surface area projected on the horizontal xy- 
plane (also denoted nominal surface area Aq). Thus in 
the Wenzel model assumption (2) takes the form 


nsY = nsL + 7lvcos6»o 


and combining this with (1) gives 

COS0O = rcos9 (3) 
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FIG. 1. Fluid in contact with a substrate with cosines corru¬ 
gation. The fluid occupy half the region where the absolute 
value of the slope is below —ta.nO. 

SO that ii 9 > it 12^ 9o > 9. 

Note that the equations above for 75 ^ and 7 gy depend 
on the assumption that the surface energies 7 sv and 7 sl 
are independent on the surface slope which may be the 
case for amorphous solids but in general not for crys¬ 
talline solids. 

Let P{s) be the probability distribution of the absolute 
value of the surface slopes s = |V/i(x)|. If the fluid make 
complete contact with the substrate thenji^l 


fluid will make complete contact with the solid (Wenzel 
state). However, if — tan0 < qho we have the situation 
shown in Fig. [T] In this case the fluid will make contact 
with the solid in half of the surface region where the 
slope s < —tan 0 and we can write 

7sl = ds P{s) (7SL + 7sv) (l + s^) 

/ OO 

ds P{s) ( 7SV (1 + s^) ^ + 7 lv ) (7) 

-tan0 ' 

and 

poo 

7sv = 7SV / ds P{s) (1 -b (8) 

Jo 

Combining (2), (7) and ( 8 ) and using (1) gives 

2 p — tanO 

cos9o “ 2 / ^ 

(9) 

If 9 is close to tt we have —tan0 Ri 0 and we can write (9) 
as 

p — tan9 

cos9o = —1 + - (cos 6 * + 1) / ds P{s) (10) 
2 Jo 


ds P{s){l + s^y^^ ( 4 ) 

For a randomly rough surface it has been shown that 

P(s) = (5) 

So 

where sq is the root-mean-square slope. For a cosines 
profile h{x,y) = hocos{qx) (where q = 27r/A, where A is 
the wavelength) one obtain 

PM = ^{2»S(6) 

for s < sov^2 and P{s) = 0 for s > so-\/2, where sq = 
h^qf ^2 is the root-mean-square slope. 

When the contact angle 9 > 7r/2, for surface roughness 
with large enough slopes, a fluid droplet may be in a 
state where the vapor phase occur in some regions at the 
nominal contact interface. This is denoted the Cassie 
state and for this case it is much harder to determine 
the macroscopic (apparent) contact angle 9q. In fact, it 
is likely that many metastable Cassie states can form. 
In general, if several (metastable) states occur, the state 
with the smallest contact angle will be the (stable) state 
with the lowest free energy. 

Let us first focus on the case where the fluid makes 
contact with a cosines profile h{x,y) = hocos(gx). In 
this case the contact will be as indicated in Fig. [TJ If 
—tan 0 is larger than the maximum slope qho = 



This equation was derived for a surface with a single 
cosines corrugation, but should hold approximately also 
for a randomly rough surface with roughness on a single 
length scale, i.e., a surface generated by the superposition 
of cosines waves (or plane waves) with equal wavelength 
but different propagation directions (in the xy-plane) and 
with different (random) phases. Such a surface will have 
the distribution of slopes given by (5). With P{s) given 
by (5) from (10) we get 


1 / 


/tan 0 A ^ 

COS 0 O = —1 H— (cos 0 -f 1 ) 

1 — exp 

2 \ 


[v So y J 


( 11 ) 


If we write 9o = tt — (fo and 9 = tt — (f with c/jq < < 1 and 
(j) « 1 we get from ( 11 ) 



Consider a randomly rough surface with roughness on 
many length scales. Most solids have surface rough¬ 
ness which is approximately self-affine fractal over several 
decades in length scales |45l|. A self affine fractal surface 
has a power spectrum which depends on the wavevec- 
tor g as a power law C{q) ~ ^}^g Hurst 

exponent 0 < iL < 1. In addition, most surfaces have 
a roll-off region for q < q^. The solid line in Fig. [2] 
shows the power spectrum of such a surface. Here we 
consider instead a rough surface composed of a discrete 
set of wavelength components with the power spectrum 
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FIG. 2. The solid line shows the power spectrum as a func¬ 
tion of the wavevector (log-log scale) of a self-afRne fractal 
surface with roughness over several decades in length scale. 
The power spectrum given by the vertical arrows (Dirac delta 
functions) correspond to a rough surface composed of a dis¬ 
crete set of wavelength components. We can consider this as 
a discretize version of the continuous power spectrum. We 
assume the separation between two nearby delta functions to 
be of order 1 decade in length scale. 



FIG. 3. The interface between a fluid droplet and a rough sub¬ 
strate. At large length scales the liquid droplet will not make 
full contact with the rough surface (Gassie state). However, 
when an apparent contact area is magnified, it is observed 
that at shorter length scales the fluid is in complete contact 
with the roughness profile. 


given by the vertical arrows (Dirac delta functions) in 
Fig. [2] We can consider this as a discretize version of 
the continuous power spectrum in Fig. [21 We assume the 
separation between two nearby delta functions to be of 
order I decade in length scale. This large separation in 
length scales allow us to integrate out, or eliminate the 
roughness which occur at length scales shorter than the 
length scale under consideration. 

We will now consider how the contact angle Oq changes 
as we gradually add more roughness components to the 
surface profile. We will use a Renormalization Group 
type of picture and study how the effective interfacial 
energy 7 sl(?t-) (and 'ysv(n)) change as we include more 
and more of the surface roughness components, i.e., we 
gradually increase n from n = 1 to the final value n = ni 
where all the roughness components are included (see 
Fig. [2]). Consider first a surface where we only include 
the shortest wavelength roughness components (associ¬ 
ated the n = 1 Dirac delta function in Fig. IlD. Weassume 
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FIG. 4. The contact angle 9o = 0(n) as a function of the 
parameter r — Atot/Ao. (a) The red symbols are the results 
obtained from the equation cos0(n) = r„cos^(n — 1) with 
n < nw (in the present case (n < nw = 7) and the blue 
symbols using (j){n + 1) = {n)/{^y2sn) for n > nw- (b) The 
limiting case where roughness occur on infinite many length 
scales (upper curve) or on a single length scale (lower curve). 


that the slope of the surface is everywhere below —tan0. 
In this case fluid will make complete contact with the 
rough surface. If we now consider the system at a lower 
magnification the surface appear smooth and flat. The 
fluid contact angle at this magnification can be calculated 
from (3) with the modified (effective) interfacial energy: 

7sl(?^ = 1) = 7SL X ri 

where ri = A(n = l)/^o is the ratio between the surface 
area A(n = 1) and the nominal surface area Aq. Using 
this effective interfacial energy, and a similar expression 
for 7sv(’T- = 1), from (3) one obtain the fluid droplet 
contact angle cos0(l) = ricos0. Let us now add the next 
shortest wavelength roughness components, represented 
by the Dirac delta function peek n = 2 in Fig. We 
assume that the slope of the surface is everywhere below 
—tan0(l). In this case fluid will make complete contact 
with the rough surface. If we now consider the system 
at a lower magnification the surface appear smooth and 
flat. The fluid contact angle at this magnification can 
be calculated from (3) with the modified (or effective) 
interfacial energy: 

7SL(n = 2) = 7sl(1) X r2 

where r 2 = A{n = 2)/A(n = 1) where A{n = 2) is the 
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surface area obtained with the n = 1 and n = 2 roughness 
components. Using this effective interfacial energy, and 
a similar expression for 7sv(’^ = 2), from (3) one obtain 
the fluid droplet contact angle cos0(2) = r2cos0(l). More 
generally, we have 

cos9{n) = r„cos0(n — 1) (13) 

where r„ = A(n)/A{n — 1). Eq. (13) is valid also for 
n = 1 if we define 9{n = 0) = 0 as the contact angle on 
the flat perfectly smooth surface, and A(n = 0) = Aq. 
If the root-mean-square slope associated with the rough¬ 
ness in each of the Dirac delta functions peeks is small 
enough, the procedure described above can be contin¬ 
ued until 9{n) is close to tt (or 180°) and we conclude 
that for a surface with roughness on many length scales 
the Wenzel state will prevail until the root-mean-square 
roughness is so high that the Wenzel equation predict a 
contact angle close to tt. We define n = nw as the small¬ 
est n for which r„.^,^,+icos0(nw) < ~1, at which point 
there is no solution to the Wenzel equation (Eq. (3)): 
cos0(nw + 1) = ?'„w-i-icos0(nw) 

Consider now adding the n = nw +1 roughness compo¬ 
nent. Since the maximum surface slope associated with 
the n = nw + 1 roughness component is larger than 
—tan0(nw), the liquid droplet will not make full con¬ 
tact with the rough surface. However, at shorter length 
scales the fluid is in complete contact with the roughness 
profile as indicated in Fig. [3l To study the development 
of the contact angle as we add more roughness we assume 
that (9) is valid, which for 9 = -k — <j) k, tt takes the form 
(12) which we now write as 

</>(« +1) = ^ (14) 

for n > nw, with (j){nw) determined by the contact angle 
in the Wenzel state which prevail when n = nw- In (14) 
Sn is the rms slope associated with the surface rough¬ 
ness contained in the n’th Dirac delta function (see Fig. 
El). Note that for a randomly rough surface s„ can be 
calculated from r„. 

We now present numerical results to illustrate the dis¬ 
cussion above. Assume that the fluid contact area on the 
flat smooth surface equals 9 = 120° as would be typical 
for Teflon. Assume that ri = r 2 = ■■ = 1.1. This corre¬ 
spond to Si = S 2 = ... ~ 0.46. Note that the total surface 
area after adding the first n roughness wavelength com¬ 
ponents is Atot/Ao = n X r 2 X ... x r„. In Fig. HKa) we 
show the contact angle 9q = 9{n) as a function of the 
parameter r = Atot/Ao. 

It is instructive to consider two limiting cases, namely 
the case where surface roughness occur over infinite 
decades in length scales (this limit can of course not be 
realized in reality as it is meaningless to consider rough¬ 
ness at length scales below the atomic dimension) and 
when roughness occur on a single length scale. In the for¬ 
mer case, we assume as above a surface roughness power 


spectrum consisting of Dirac delta function peaks sepa¬ 
rated by ~ 1 decade in lengthscale. We consider a surface 
with a finite rms slope and we assume that each Dirac 
delta function contribute equally to the rms slope. Hence 
since there are infinite many Dirac delta functions each 
of them must have an infinitesimal weight corresponding 
to roughness with infinitesimal amplitude. It is clear that 
in this case the condition r„.^,^,_|_icos 0 (nw) < — 1 will be 
satisfied only when the contact angle 0 (nw) differ from 
TT by an infinitesimal amount, i.e., the Wenzel state will 
prevail (minimum free energy state) until r is so large 
that Eq. (3) predict 9 = tt after which the Cassie state 
will prevail (with 9 = tt). 

For the second limiting case of roughness on a single 
length scale we use ( 11 ) to estimate the contact angle 
6*0 for the Cassie state. In Fig|4Kb) we show the contact 
angles for these two limiting cases. 

We now present results for the contact angle as a func¬ 
tion of the roughness parameter r based on an exact 
numerical treatment of an one-dimensional (ID) model. 
Using standard procedures (see, e.g., Appendix D in Ref. 
[ 4 ^) we have generated randomly rough ID surfaces and 
studied the fluid contact angle using the method de¬ 
scribed in Ref [s^. We consider surfaces with rough¬ 
ness extending over what we denote as a “narrow” and 
a “wide” range of length scales. The roughness power 
spectral density (PSD) C{q) of the two types of surfaces 
are shown in Fig. [S] Both surfaces have the rms slope 
0.25 and Hurst exponent H = 0.8, but different large 
wavevector cut-off. As a result, the rms roughness is 
larger for the surface with the more narrow PSD. This is 
illustrated in Fig. [6] which shows the height distribution 
Ph of both surfaces. 

We will show below that in spite of the larger rms 
roughness amplitude of the surface with the more nar¬ 
row PSD, the (apparent) contact angle is largest on the 
surface with the widest PSD. This shows that the rms 
roughness is irrelevant for the contact angle which in¬ 
stead is determined by the rms slope and how different 
length scales contribute to the rms slope. 

Fig. [7] shows the contact angle 6*0 as a function of the 
parameter r = Atot/Ao obtained from numerical simula¬ 
tions. Note that the surface with the more wide distri¬ 
bution of roughness lengthscales has larger contact angle 
in spite of the fact that both surfaces have the same rms 
slope. In general, several metastable droplet states are 
possible, and which state occur depends on the prepa¬ 
ration procedure. We note, however, that the minimum 
free energy state is the state with the lowest contact an¬ 
gle. Thus, the Cassie states observed for r < 2 in Fig. 
[7] are metastable states, and the ground state is in this 
case the Wenzel state. In the theory presented above we 
assumed that the system was in the minimum free energy 
state and no Cassie state could occur above the Wenzel 
line for r < 2 . 

The result in Fig. [7] is consistent with the theory pre- 
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FIG. 5. The surface roughness power spectra C{q) of surfaces 
with the rms slope 0.25 and Hurst exponent H = 0.8 but 
different large wavevector cut-off. 



h (gm) 


FIG. 6. The height distribution Ph of surfaces with the rms 
slope 0.25 and Hurst exponent H — 0.8 but different large 
wavevector cut-off. 


diction presented in Fig. 21 In the numerical simulations 
the contact angle in the Cassie state is smaller than in the 
theory, but this reflects the fact that the roughness occur 
over a wider range of length scales in the theory as com¬ 
pared to the numerical model. This is indeed confirmed 
by numerical calculations presented by the authors in 



r=A„ /A„ 

FIG. 7. The contact angle 6o as a function of the parameter 
r = Atot/Ao obtained from numerical simulations. The solid 
line is given by the Wenzel equation (3). 


(see Fig. 10 therein) where the surface roughness was 
characterized by a much larger number of length scales 
covering about 3 decades. 

To summarize, we have shown when the number of 
length scales of roughness increases, the transition from 
the Wenzel state to the Cassie state approach the thresh¬ 
old value r = ^tot/Tlo = — l/cos0, and the contact angle 
in the Cassie state approaches 180°. 
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